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Abstract 



Dynamic heterogeneity is now recognised as a central aspect of structural relaxation in 
disordered materials with slow dynamics, and was the focus of intense research in the 
last decade. Here we describe how initial, indirect observations of dynamic heterogene- 
ity have recently evolved into well-defined, quantitative, statistical characterisations, 
in particular through the use of high-order correlation and response functions. We 
highlight both recent progress and open questions about the characterisation of dy- 
namic heterogeneity in glassy materials. We also discuss the limits of available tools 
and describe a few candidates for future research in order to gain deeper understanding 
of the origin and nature of glassiness in disordered systems. 



Fig. 0.1 Three examples of dynamical heterogeneity. In all cases, the figures high- 
light the clustering of particles with similar mobility. (Left) Granular fluid of ball 
bearings, with a colour scale showing a range of mobility increasing from blue to 
red ( |Keys et al., 2007| . (Centre) Colloidal hard sphere suspension, with most mobile particles 
highlighted ( | Weeks et al., 2000[ ). (Right) Computer simulation of a two-dimensional system 
of repulsive disks. The colour scheme indicates the presence of particles for which motion is 
reproducibly immobile or mobile, respectively from blue to red (Widmcr-Cooper et al., 2008). 

0.1 Introduction 

0.1.1 Dynamical heterogeneity in glassy materials 

The glass transition is often cited as a profound outstanding problem in condensed 
matter physics. This field may be contrasted with that of simple liquids, for which the 
broad picture is now well-established, and appropriate theoretical methods are well- 
developed ( jHansen and McDonald, 1986[ ) . Why, then, is the glass problem so difficult? 

From a theoretical perspective, a central difficulty arises from the importance of 
fluctuations in glassy systems. Both the liquid and the glass have disordered structure, 
so even if all molecules in the system are identical, they experience different local 
environments. In the liquid, these differences can be neglected: one may infer the 
behaviour of the system from that of a typical particle in a typical environment. Thus, 
for example, microscopic properties, such as the rate with which particles diffuse in 
the liquid, arc directly related to bulk properties, such as the viscosity. However, as the 
glass transition is approached, it becomes increasingly difficult to characterise 'typical' 
particles and 'typical' environments because a variety of different behaviors emerges. 
Within a given interval of time, some particles may move distances comparable to their 
size, while others remain localised near their original positions. Thus, on these time 
scales, we can refer to 'mobile' and 'immobile' particles. Of course, on long enough 
time scales, ergodicity ensures that particles become statistically identical. 

In case this discussion seems rather abstract, we refer the reader to Fig. 10.11 Here, 
we show several systems in which mobile and immobile particles can be identified in 
particular trajectories using different methods. Strikingly, these images reveal that par- 
ticles with different mobilities do not appear randomly in space but are clustered. This 
observation suggests that structural relaxation in disordered systems is a nontrivial 
dynamical process. In its narrow sense, the term 'dynamical heterogeneity' encapsu- 
lates the spatial correlations observed in Fig. 10.11 However, the term is frequently 
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used in a broader sense, referring to a range of fluctuation phenomena that arise from 
deviations from the 'typical' behaviour ( |Ediger, 2000[ ). 

Over the last decade, it has become clear from experiments and computer simula- 
tions that a variety of glassy systems display the kind of clusters shown in Fig. 10.11 
Experimentally, their existence can be inferred from experiments in molecular liq- 
uids dEdiger, 2000l|Sillescu, 1999l|Richert, 2002[|R7chert et al, 2010| ), while direct ob- 
servation of single particle motion makes them vivid in colloids and granular me- 
dia ( [Weeks et al, 2000|Keys et al, 2007|Kegel and van Blaaderen, 2000|Dauchot et al, 2005 1 



In computer simulations, spherical particles with simple pair potentials have been used 
as models for both colloidal and molecular systems, with dynamically heterogeneous 
behaviour clearly present in a variety of models. Dynamical heterogeneity has also 
been investigated in a large number of more schematic lattice models, such as kinct- 



ically constrained or lattice glass models (Ritort and Sollich, 20031. This important 



body of experimental and computational observations has also stimulated important 
theoretical developments since they represent a new set of observations against which 
existing theories can be confronted. 

So far, the reader may be unconvinced of the difficulty of the problem. After all, 
theoretical methods for analysing systems with large fluctuations and correlated do- 
mains already exist: the methods developed to describe critical phenomena, such as 
field theory and the renormalisation group. These ideas are indeed central to this 
chapter, but their application to glassy liquids has required substantial new insight 
into the nature of the relevant fluctuations and observables. The reason is that the 
distinction between mobile and immobile particles is in essence dynamical. Therefore, 
if one analyses static snapshots of viscous liquids, there is little evidence of increasing 
fluctuations or heterogeneity, at least when analysed using standard liquid state corre- 
lation functions. Instead of an ensemble of snapshots, one must apply the methods of 
critical phenomena to an ensemble of 'movies' (i.e. trajectories, or dynamical histories) 
of the system. This will be the approach that we will follow in later sections. 



0.1.2 Suitable probes for the emergence of glassiness 

Returning to Fig. 10.11 the observation of clusters of mobile particles (or at least of 
regions with correlated mobility) raises many questions. What is the nature of the 
clusters? Indeed, are they of the same nature in each case? What determines their size? 
How is their size distribution related to the relaxation time of the system, if at all? 
Do the particles in fast or slow clusters have different local environments that can be 
characterised by any simple structural measure? From a theoretical point of view, it is 
natural to ask whether the correlated regions represent the cooperatively rearranging 
regions imagined long ago by Adam-Gibbs ( |Adam and Gibbs, 1965D ; whether they 
mirror some soft elastic modes in the system; how they might connect to locally- 
ordered domains; or whether they reveal the presence of localised defects that facilitate 
structural rearrangement. 

However, before addressing these ambitious questions, we must answer a more 
prosaic one. How can the 'clusters of mobile and immobile particles' be defined and 
characterised? In Fig. 10. 1[ the mobile particles were identified by a threshold on their 
displacement, over a particular time scale. Indeed, many early papers have suggested 
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different ways to define 'mobility' and 'clusters' that are similar in spirit but different 
in practice, as we detail in the following sections. However, one must certainly evaluate 
how strongly the cluster properties depend on thresholds or time scales, which might 
hinder firm conclusions and may prevent fair comparison between different systems. 

In more recent years, it has become possible to define and measure observables that 
can be determined without arbitrariness in a range of systems, are amenable to analytic 
theory and scaling arguments, and may sometimes be inferred from experimental data 
even in molecular liquids. These are known as 'four-point correlation functions' and 
are now broadly accepted as standard tools for analysing dynamical heterogeneity. 

Within this toolbox, a central role is played by the four-point dynamical susceptibil- 
ity X4(i). Loosely speaking, it measures the number of particles involved in correlated 
motion on times scales of the order of t. To interpret the dynamical susceptibility, it 
is useful to invoke an analogy with critical phenomena, in which a (static) correla- 
tion length £ diverges at a critical temperature T c , accompanied by the spontaneous 
appearance of an order parameter. This divergence is associated with a diverging sus- 
ceptibility Xj which may be measured either through the fluctuations of the order 
parameter or through the response of the order parameter to its conjugate field. When 
considering such phenomena in glasses, a problem arises, in that a static order pa- 
rameter and its conjugate field are not known. Instead, a fairly good dynamical order 
parameter is given by any generic dynamical two-point correlator, e.g. density-density, 
displaying the slowing down of the dynamics. This motivated the use of a four-point 
dynamical susceptibility associated to spontaneous fluctuations of the dynamical or- 
der parameter. This allows the identification of a dynamic length scale, ^(t), and a 
susceptibility, Xi(t)i by analogy with conventional critical phenomena. 

Returning to Fig. I0.I[ simulation studies and recent experiments indicate that the 
clustering of mobile particles is directly linked with an increasing susceptibility Xi 
and an increasing correlation length scale £4. The central part of this chapter will 
be devoted to a discussion of these four-point functions. We will discuss some of the 
insights that they have revealed into the nature of glassy behaviour in liquids, col- 
loids, and granular media. However, these studies also revealed that interpretation of 
four-point functions may be somewhat ambiguous, while direct measurements of corre- 
lation length scales remain difficult. Additionally, the averaging procedure inherent in 
the four-point functions means that they may obscure important features of dynamic 
heterogeneity such as the cluster shape and the nature of interfaces between clusters of 
mobile and immobile particles. Towards the end of the chapter, we will discuss a range 
of alternative observables that complement the information available from four-point 
functions. 



0.2 Observables for characterising dynamical heterogeneity 

0.2.1 Two-point observables and their inadequacy 

We begin with a review of some two-point functions that are used to characterise 
simple liquids. We will show that these functions are largely blind to the dynami- 
cally heterogeneous behaviour shown in Fig. 10. 1\ motivating the discussion of more 
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discriminating observable^ 

For any fluid of particles, a natural quantity to measure is the structure factor, 



S(q) = ^ (PqWP-qW) , (0.1) 



where brackets indicate an ensemble average, and the Fourier component of the density 
is 

N 

PqW=E e4q ' ri(t) ' ( - 2 ) 

1=1 

with N being the number of particles and r j (t) being the position of particle i at time 
t. The structure factor gives information about the strength of density fluctuations 
on a length scale 27r/|q|. However, its behaviour in the vicinity of the glass transition 
is unremarkable, with no hint of the dynamic clustering of Fig. 10.11 Although static 
heterogeneities in the density would directly imply the existence of dynamic hetero- 
geneity, the reverse is not true. Thus, dynamic heterogeneity related to the motion of 
particles has a much more subtle origin. Note that more complicated static correlation 
functions have been studied ( |Debenedetti, 1996] ), especially in numerical work, all at- 
tempting to identify fluctuations of some prescribed sort of local order (translational, 
orientational, etc.). Until now, there are no strong indications of a diverging, or even 
substantially growing, lengthscale ( Menon and Nagel, 1995 , Fer nandez et ai, 2006 ). In 



order to sec the growth of a static amorphous correlation, one possibility is to intro- 
duce the so-called "point-to-set" correlation function that we discuss later, see section 

We therefore turn to dynamical observables. A quantity relevant for light and 
neutron scattering experiments is the intermediate scattering function, 



F(q,t) = l<p q (*)p_ q (0)>. (0.3) 



Measurements of this function by neutron scattering in supercooled glycerol flWuttke et ai, 1996[ ) 
are shown for different temperatures in Fig. 10.21 These curves suggest a first, rather 
fast, relaxation to a plateau followed by a second, much slower, relaxation. The plateau 
is due to the fraction of density fluctuations that are frozen on intermediate timescales, 
but eventually relax during the second relaxation. The latter is called 'alpha-relaxation', 
and corresponds to the structural relaxation of the liquid. The plateau is akin to 
the Edwards- Anderson order parameter, defined for spin glasses which measures the 
fraction of frozen spin fluctuations ( |Binder and Kob, 2005[ ). Note that the Edwards- 
Anderson parameter continuously increases from zero below the critical temperature 



in the conventional spin glass transition (Mczard et ai, 19881, while for structural 



glasses, a finite plateau value seems to emerge above any putative transition. 

The full decay of the intermediate scattering function can be measured only within 
a relatively small range of temperatures. In order to track the dynamic slowing down 

1 There are other examples of systems for which standard two-point correlators are blind to inter- 
esting intermittent (or heterogeneous) effects, such as turbulence or financial markets, and for which 
higher order correlations are informative. 
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Fig. 0.2 Temperature evolution of the normalised intermediate scattering function, 
4>q(i) = S(q, t)/S(q, 0), for supercooled glycerol ( |Wuttke et al, 1996[ ). Temperatures de- 
crease from 413 K to 270 K from left to right. The lines are fits with a stretched exponential 
form. 

from microscopic to macroscopic timescales, other correlators have been studied. A 
very popular one is the measurement of the dielectric linear susceptibility which can 



be followed over up to 18 decades of frequency ( |Lunkenheimer and Loidl, 200 2). It is 
generally accepted that different dynamic probes reveal similar temperature depen- 
dences for the relaxation time, at least as long as the probes measure local motion. 
In broad terms, the essential features on supercooling are a dramatic increase in the 
correlation time, and a broad distribution of time scales in the system characterised 
in the time-domain by non-cxponcntial relaxation functions. 

It is increasingly accepted that the presence of a such broad distributions of time 
scales in glassy systems is associated with the presence of mobile and immobile do- 
mains. However, the size and shape of these domains, or even their very existence, 
cannot be deduced directly from F(q,t). We are therefore motivated to consider more 
advanced correlation functions. 



0.2.2 Indirect evidence: Intermittency and decoupling phenomena 

We now return to the presence of mobile and immobile particles in the supercooled 
phase. Simulation studies are ideal for studying fluctuation properties, since accurate 
trajectories for all particles are accessible, over long time scales. These features are dis- 
played in Fig. 10.31 which shows that while the averaged mean-squared displacements 
are smooth functions of time, time signals for individual particles clearly exhibit spe- 
cific features that are not observed unless dynamics is resolved both in space and time. 
In this figure, we observe that particle trajectories are very intermittent, being com- 
posed of a succession of long periods of time where particles simply vibrate around 
well-defined locations, separated by rapid 'jumps'. Vibrations were previously inferred 
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Fig. 0.3 Time resolved squared displacements of individual particles in a simple model of a 
glass-forming liquid composed of Lennard- Jones particles (Bcrthicr and Kob, 2007). The av- 
erage is shown as a smooth full line and time is expressed in units of structural relaxation time 
T a . Trajectories are composed of long periods of time during which particles vibrate around 
well-defined positions, separated by rapid jumps that are statistically widely distributed in 
time, underlying the importance of dynamic fluctuations. 

from the plateaux observed at intermediate times in the mean squared displacements 
or intermediate scattering functions, but the existence of jumps that are statistically 
widely distributed in time cannot be revealed from averaged quantities only. The fluc- 
tuations in Fig. 10.31 suggest, and direct measurements confirm, the importance played 
by fluctuations around the averaged dynamical behaviour to understand structural 
relaxation in glassy materials. 

Remaining at the single particle level, these fluctuations can be characterised 
through the (time-dependent) distribution of particle displacements. This is the self- 
part of the van-Hove function, defined as 



G * ( r ' *) = ( E ^ - ^ (*) - * (0)] ) ) ■ (0-4) 



For an isotropic Gaussian diffusive process, one has G s (r, i) ~ cxp (— t \ . While 

simple liquids are well-described by such a distribution, simulations of glassy systems 
instead reveal strong deviations from Gaussian behaviour on the timescales relevant for 



structural relaxation (Kob et at, 1997). In particular they reveal 'broad' tails in the 
distributions that are much wider than expected from the Gaussian approximation. 
These tails are in fact well described by an exponential, rather than Gaussian, decay 
in a wide time window comprising the structural relaxation, 

G s (r,t)~exp(-Jjl) , (0.5) 
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Fig. 0.4 Decoupling between viscosity (full line) and self-diffusion coefficient (symbols) in 
supercooled ortho-terphenyl ( |Mapes et al, 2006[ ). The dashed line shows a fit with a 'frac- 
tional' Stokes-Einstein relation, D s ~ (T/rf) with £ ~ 0.82 instead of the expected value 

C = i- 

which is a direct consequence of the intermittent motion shown in Fig. l0.3l (Chaudh uri et al., 2007[ ) 
These tails reflect the existence of a population of particles that moves distinctively 
further than the rest and appears therefore to be much more mobile. This observation 
implies that relaxation in a viscous liquid differs qualitatively from that of a normal 
liquid where diffusion is close to Gaussian, and that a non-trivial statistics of single 
particle displacements exists in materials with glassy dynamics. 

Another influential phenomenon that was related early on to the existence of dy- 
namic heterogeneity is the decoupling of self-diffusion (D s ) and viscosity (77). In the 
high temperature liquid, self-diffusion and viscosity are related by the Stokes-Einstein 
relation ( Hansen and McDonald, 1986[ ), D s rj/T = const. For a large particle moving 



in a fluid the constant is equal to l/(6nR) where R is the particle radius. Physi- 
cally, the Stokes-Einstein relation means that two different measures of the relaxation 
time, R 2 /D s and rjR 3 /T, lead to the same timescale up to a constant factor. In su- 
percooled liquids this phenomcnological law breaks down, as shown in Fig. 10.41 for 
ortho-terphenyl ( |Mapes et al., 2006 ). It is commonly found that D' 1 does not increase 



as fast as r\ so that, at T g , the product D s r\ has significantly increased as compared 
to its Stokes-Einstein value. The Stokes-Einstein 'violation' factor is larger for more 
fragile liquids, and can be as high as 10 3 . This phenomenon, although less spectacu- 
lar than the overall change of viscosity, is a strong indication that different ways to 
measure relaxation times lead to different answers and, thus, is a strong hint of the ex- 



istence of broad distributions of relaxation timescales (Stillinger and Hodgdon, 1994 



Tarjus and Kivelson, 1995) 



Indeed, a natural explanation of this effect is that different observables probe the 
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underlying distribution of relaxation times in different ways ( |Ediger, 2000[ ) . For exam- 
ple, the self-diffusion coefficient of tracer particles is dominated by the more mobile 
particles whereas the viscosity or other measures of structural relaxation probe the 
timescale needed for every particle to move. An unrealistic but instructive example is 
a model where there is a small, non-pcrcolative subset of particles that are blocked 
forever, coexisting with a majority of mobile particles. In this case, the structure 
never fully relaxes but the self-diffusion coefficient is non-zero because of the mo- 
bile particles. Although unrealistic since all particles move in a viscous liquid, this 
example shows how different observables are likely to probe different moments of 
the distribution of timescales, as explicitly shown within several theoretical frame- 



works flTarjus and Kivelson, I995l|Jung et ai, 200 4, Hedg es et al., 2007l|Heuer, 2008| >. 



0.2.3 Early studies of dynamic heterogeneity 

The phenomena described above, although certainly an indication of spatio-temporal 
fluctuations, do not allow one to study how these fluctuations are correlated in space. 
However, this is a fundamental issue both from the experimental and theoretical points 
of view, as discussed in the introduction. To discriminate between different explana- 
tions of glassy behaviour, it would be useful to know: How large are the regions that 
are faster or slower than the average? How does their size depend on temperature? 
Arc these regions compact or fractal? 

These important questions were first addressed in pioneering works using four- 
dimensional NMR ( |Tracht et ai, 1 99 8,Roins berg et ai, 2001[ ), and by directly probing 
fluctuations at the nanoscopic scale using microscopy techniques. In particular, Vidal 
Russcl and Isracloff used Atomic Force Microscopy techniques ( |Vidal Russell an"d~I sraeloff, 2000[ ) 
to measure the polarisation fluctuations in a volume of size of few tens of nanometers in 
a supercooled polymeric liquid (PVAc) close to T g . In this spatially resolved measure- 
ment, the hope is to probe a small enough number of dynamically correlated regions, 
and to detect their dynamics. Indeed, the time signals shown in Ref. flVidal Russell and Israeloff, 2000[ ) 
show a very intermittent dynamics, switching between moments with intense activity, 
and moments with no dynamics at all, suggesting that extended regions of space in- 
deed transiently behave as fast and slow regions. A much smoother signal would have 
been measured if dynamically correlated 'domains' were not present. 

Spatially resolved and NMR experiments arc quite difficult. They give undisputed 
information about the typical lifetime of the dynamic heterogeneity, but their determi- 
nation of a dynamic correlation lengthscale is rather indirect, and has been performed 
on a small number of liquids in narrow temperature windows. Nevertheless, an agreed 
consensus is that these experiments reveal the existence of a non-trivial dynamic cor- 
relation length emerging at the glass transition, where it reaches a value of the order 
of 5 to 10 molecule diameters ( |Ediger, 2000[ ). 

More recently, studies of systems for which single-particle dynamical data is avail- 
able have led to direct observation of clusters of mobile particles, such as those shown 
in Fig. 10. II This is possible in materials where the particles are big enough to be 
visualised through microscopy, such as colloids, or with a camera, such as granular 
materials. In early studies, mobile and immobile particles were usually identified by 
using a threshold on their displacement within a given time interval, dividing the par- 
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Fig. 0.5 Early measurements of dynamical heterogeneity from computer simulation. 
(Left) The main figure shows the probability that a particle is within the active population, 
given that it is a distance r from another active particle ( |Kob et al., 1997] ). The increasing 
correlation between active particles is apparent on decreasing the temperature. (Right) The 
length distribution of string-like clusters of active particles ( |Donati et al., 1998| . 

tides into sub-populations. This then allows characterisation of correlations within the 
populations, as shown in Fig. 10.51 For example, ratios of radial distribution functions 
give the probability that particles in the vicinity of a mobile particle are themselves 
mobile. Other work concentrated for instance on the morphology and size of clusters 
of mobile particles, as shown in Fig. 10.51 

Results such as those of Fig. 10.51 gave clear evidence of large fluctuations and dy- 
namical heterogeneity. However, unambiguous identification of a mobile population 
of particles proved difficult, with most distributions of mobility showing broad but 
unimodal distributions. Similarly, the identification of connected clusters of mobile 
particles introduces further ambiguity into the data processing, with their size and 
shape depending quite strongly on the definition of the mobile population. To achieve 
unambiguous and system-independent definitions of length scales and correlation vol- 
umes, it is useful to define correlation functions that do not involve separating particles 
into distinct populations, nor the identification of connected clusters. Four-point func- 
tions are a natural choice in this regard, and will be discussed in the next chapter. 

0.2.4 Higher order correlations: four-point functions 

Definitions. In the previous section, we considered the probability that a particle 
is within some mobile population, given that it has a mobile particle a distance r 
away (recall Fig. 10. 5|) . While the measurement of such a probability requires the iden- 
tification of a mobile population, there is a straightforward alternative which con- 
tains similar information. We first define a continuously varying 'mobility' Cj(i,0) 
which indicates how far or how much particle i moves between times t = and 
t. Then, given two particles at separation r, one can measure the degree to which 
their mobilities are correlated. To this end, it is convenient to define a 'mobility field' 
through ( |Bennemann et al, 1999l|Donati et al, 1999l|Glotzer et al, 2000"] ) 

c(r;t,0)=$^Ci(t,0)*(r-ri). (0.6) 

i 
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Then, the spatial correlations of the mobility are naturally captured through the cor- 
relation function ( |Dasgupta et at, 1991[ ) 

G 4 (r; t) = (c(r; t, 0)c(O; t, 0)) - (c(r; t, 0)) 2 , (0.7) 

which depends only on the single time t and the single distance r = |r| as long as 
the average is taken at equilibrium in a translationally invariant system. The analogy 
with fluctuations in critical systems becomes clear in Eq. (|0.7[) if one considers the 
mobility field c(r; t, 0) as playing the role of the order parameter for the transition, 
characterised by non-trivial fluctuations and correlations near the glass transition. 

Often, the mobility Ci(t, 0) is itself a two-point function. For example, to measure 
mobility on a length scale 2ir/q, one might consider Oi(q,t) = e lq ' ri( -*' ) and Ci{t, 0) = 
Oi(q, t)oi(— q, 0). In this case, Oi(q, t) is related to a Fourier component of the density 
of the system, and the average of q (t, 0) is the self-part of the intermediate scattering 
function F(q, t) defined in Eq. (|0.3[) . Moving from a particle observable Oj(q, t) to a 
field o(r; q, t), one arrives at 

G 4 (r;t) = ( O (r;q,t) O (r;-q,0)o(0;q,t) O (0;-q J 0)) - (o(r; q, t)o(r; -q, 0)) 2 . (0.8) 

This correlation function is quartic in the operator o, so it is known as a 'four-point 
function'. It measures correlations on a length scale r, associated with motion between 
time zero and time t; it depends additionally on the length scale q -1 used in the 
definition of the particle mobility Ci(t, 0). Since structural relaxation typically involves 
particle motion over a distance comparable to the particle size R, one typically chooses 
q ~ 1/R and studies the remaining t and r dependences. 

This definition of a real-space correlation function of the mobility represents a vital 
advance in the characterisation of dynamical heterogeneity. In particular, it allows the 
language of field theory and critical phenomena to be used in studying dynamical 
fluctuations in glassy systems. By analogy with critical phenomena, if there is a single 
dominant length scale £4 then one expects that for large-r, the correlation function 
decays as 

G4(r;t)~^e-'/W*>, (0 . 9) 

with p an exponent whose value is discussed below. It is also natural to define the 
susceptibility associated with the correlation function 



Xi {t) = J drGi{r-t). (0.10) 

If the prefactor A(t) were known, the susceptibility x 4 (t) could be used to extract 
the typical number of particles involved in correlated motion. That is, x±(t) m ay be 
interpreted as the size of the correlated clusters in Fig. 10.11 

Further, Xi(t) can a l so be obtained from the fluctuations of the total mobility 
C(t,0) = / d d rc(r;t,0), through 

X4 (t) = N[(C(t,0f) (C(t,0)n (0.11) 

In practice, this formula allows an efficient measure of the degree of dynamical het- 
erogeneity, at least in computer simulations and in experiments where the dynamics 
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can be spatially and temporally resolved. As long as the observable c(r; t, 0) is chosen 
appropriately, this observable can be measured in a wide variety of systems, and serves 
as a basis for fair comparisons of the extent of dynamical heterogeneity. 

The spin glass perspective: four-point functions in space and time. It is interesting 
to note that four-point functions have their origin in spin glass physics, where they 
were used to investigate the onset of long-ranged amorphous order. The key insight 



is due to Edwards and Anderson (Edwards and Anderson, 19751. In a spin glass, one 



considers a set of N localised degrees of freedom, the spins, which we denote by s x , with 
x denoting the position in space. The spins interact by quenched random couplings. 
Then, two-point correlation functions between spins such as ^2 x (s x s x+r ) typically 
vanish for r ^ 0, since sites separated by a distance r may be either correlated or 
anti-correlated, with equal probability. The Edwards- Anderson solution is to consider 
instead the static spin-glass correlation xsg = ^ 1 Sz^Sz+r) 2 which receives a 
positive contribution for both correlated and anti-correlated sites. 
If one then considers a four-point dynamic function such as 

G4 ( r; *) = N m< S *( K+r(0)s K (i)Sx+rW}, (0.12) 
x 

then this approaches the static spin-glass correlation at long times. In spin glasses, 
spontaneous symmetry breaking occurs into an amorphous solid with long-ranged 
order at a second order spin glass transition. Near such a phase transition, G±{r;t) 
obeys critical scaling behaviour similar to that shown by static correlation functions 
near familiar phase transitions, and the corresponding susceptibility Xi{t) diverges. 

In structural glasses, the possibility of a static transition to an amorphous ordered 
state remains highly controversial. However, four-point functions such as G±{r;t) and 
Xi{t) can be usefully employed in the fluid state to characterise dynamical fluctuations, 
regardless of their possible connection to any critical point. 

Four-point susceptibilities in molecular, colloidal and granular glasses. The function 
Xi{t) has been measured in computer simulations of many different glass- forming liq- 



uids, by molecular dynamics, Brownian and Monte Carlo simulations, see e.g. (Bennemann et al., 1999 
Donati et al, 1999|Glotzer et al., 2000IFranz et al. 1999|Lacevic et al, 2003|Berthier, 2004 



Vogel and Glotzer, 2004|Berthier, 2007|Parisi, 1999|Parsaeian and Castillo, 2008| ). An 



example is shown in Fig. 10.61 for a Lennard- Jones liquid. The qualitative behaviour is 
similar in all cases ( |Franz and Parisi, 2000|Toninelli et al., 2005|Berthier et al., 2007a| ): 
as a function of time Xi{t) nrs t increases, it has a peak on a timescale that tracks the 
structural relaxation timescale and then it decreases. As mentioned above, the de- 
crease of Xi{t) at long times constitutes a major difference with spin glasses. In a spin 
glass, Xi would be a monotonically increasing function of time whose long-time limit 
coincides with the static spin glass susceptibility. Physically, the difference is that the 
spin glass state which emerges at the transition is critical or marginally stable, i.e. 
characterized by singular static responses. 

The peak value of Xiif) measures the volume over which the structural relaxation 
processes are correlated. Therefore, the most important result obtained from data such 
as presented in Fig. IO.6l is the temperature evolution of the peak height, which is found 
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Fig. 0.6 Time dependence of Xi{t) quantifying the spontaneous fluctuations of the self- in- 
termediate scattering function in a molecular dynamics simulation of a Lennard- Jones super- 
cooled liquid ( |Berthier, 2004| ). For each temperature, Xi(t) nas a maximum, which shifts to 
larger times and has a larger value when T is decreased, revealing the increasing lengthscale 
of dynamic heterogeneity in supercooled liquids approaching the glass transition. 

to increase when the temperature decreases and the dynamics slows down. From such 
data, one gets direct evidence that the approach to the glass transition is accompanied 
by the development of increasingly long-ranged spatial correlations of the dynamics. 
Note that if the dynamically correlated regions were compact, the peak of %4 would be 
proportional to £f (in d spatial dimensions), thus directly relating Xi{t) measurements 
to that of the relevant lengthscale £4 of dynamic heterogeneity 



In experiments, direct measurements of Xi{t) have been made in colloidal ( Weeks et gl, 2007 ) 
and granular materials ( |Dauchot et gl., 2005 , Ke ys et gl, 2007 ) close to the colloidal 
and granular glass transitions, and also in foams QMayer et a/., 2004] ) and gels flDuri and Cipelletti, 2006[ ), 
because dynamics is more easily spatially and temporally resolved in those cases. The 
results obtained in all these cases are again broadly similar to those shown in Fig. 10.61 
both for the time dependence of X4(£) and its evolution upon a change of the relevant 
variable controlling the dynamics. 

A major issue is that obtaining information on the behaviour of Xi{t) an d Gi(r; t) 
from experiments on molecular systems is difficult. In molecular liquids, it remains a 
difficult task to resolve temporally the dynamics at the nanometer scale. Such mea- 
surements are however important because numerical simulations and experiments on 
colloidal and granular systems can typically only be performed for relaxation times 
spanning at most 5-6 decades. On the other hand, in molecular liquids, up to 14 decades 
can be measured, and extrapolation of simulation data all the way to the experimental 
glass transition is fraught with difficulty. Indirect estimates of Xi(t) from experiments 
will be discussed below. 

These results are also relevant because many theories of the glass transition assume 
or predict, in one way or another, that the dynamics slows down because there are 
increasingly large regions over which particles have to relax in a correlated or cooper- 
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ative way, see Sec. 10.3.21 and Ref. flTarjus, 2010D . However, this lengthscale remained 
elusive for a long time. Measures of the spatial extent of dynamic heterogeneity, in 
particular Xi{t) an d Gi{r;t), seem to provide the long-sought evidence of this phe- 
nomenon. This in turn suggests that the glass transition can indeed be considered as 
a form of critical phenomenon involving growing time scales and length scales. This is 
an important progress towards the understanding of the glass transition phenomenon, 
even though a clear and conclusive understanding of the relationship between dynamic 
lengthscalcs and relaxation timescales is still the focus of intense research. 

Dependence of x<i(t) on the observable and probe length scale. As discussed above, 
one may define a four-point function Gi(r;t) starting from any suitable mobility 
c(r;i, 0). Indeed, many candidates have been considered. It is not our intention to 
give a detailed list, but a few comments are in order. 

A natural choice for Xi{t) 1S to start from Eq. (|0. 1 1|) and to take C(t, 0) = 
p q (t)p_q(0) as the autocorrelation of a single Fourier component of the density. In this 
case, the average of C(t, 0) is the intermediate scattering function F(q, t) of Eq. (|0.3[) . 
In computational studies, it is often more convenient to construct instead C(t, 0) from 
a simple sum over particles. That is, one defines the single-particle mobility 

/ i (q )i ,0) = e ^C r *W- r ^°)) ! (0.13) 

whose average is the self-part of F(q,t). The real-space four-point function is then 
given by Eq. (|0.7[) . and the definition of Xi{t) follows. These two definitions of Xi{t) are 
not equivalent. Differences between them were discussed in Ref. ( |Lacevic et al, 2003] ), 
where it was concluded that they contain similar information. 

Physically speaking, the key point is that as particle i moves away from its initial 
position r.j(0), the function /»(q; t, 0) decays from a value of unity, reaching zero when 
the particle has moved a distance of order (7r/|q|). Once the particle has moved fur- 
ther than this, the oscillations in the cosine function imply that averages of fi give 
numbers close to zero. Based on this physical interpretation, other choices for a(t, 0), 
including step functions, or smoothly decaying functions were used ( |Franz et al, 19991 
|Bennemann et al, 1999|[Lacevic et al, 20 03 } Bcrthier , 2004[ ). As expected on physical 
grounds, constructing four-point functions based on these choices for Cj(t, 0) again 
leads to qualitatively similar behaviours. 

Yet another choice is to use a function Ci(t, 0) that depends not just on the positions 
at time zero and time t, but also on the whole history of the particle between these 
times. In particular, one may take a 'persistence' function which takes a value of 
unity if the particle remains within a distance a of its initial position for all times 
between and t; otherwise it takes the value zero. Again, one observes a broadly 
similar behaviour ( |Chandler et al, 2006[ ). 

All these choices for the definition of the local mobility Cj(i,0) involve a 'probe 
length scale', which is fixed by the choice of measurement, in contrast to the sought 
dynamic length scale £/t(t) which should be a physical property of the system. While 
the specific form of c(r; t, 0) is typically unimportant for the qualitative behaviour of 
four-point functions, the quantitative results do depend on the probe length scale. 

Typically, if the probe length scale a is of the order of the particle diameter or 
smaller, c(r; t, 0) measures local motion, and this is often the scale on which heterogene- 
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ity is most apparent. As the probe length scale is increased, contributions to Xi(t) come 
from pairs of particles that remain correlated over distances comparable to a and, typ- 
ically, such correlations weaken as a increases, reducing x±{t) {Chandle r et al, 2006[ ). 
Similarly, reducing a also reduces Xi(t) as short-scale motion corresponding to ther- 
mal vibrations are also typically uncorrclated. Therefore, Xi(t) is usually maximal 
for a probe length scale comparable to the particle size, and it is fixed to a constant 
when comparing data at different temperatures or densities. An alternative choice is 
to adjust the probe lengthscale a at different state points such that Xi{t\ a ) reaches its 
absolute maximum, this can be very important for some systems like granular media 
close to the rigidity transition where the maximum is reached for values of the probe 



length far below the particle size ( |Lechenault et al, 2008a||Lechenault et al., 20086 
|Heussinger et al., 2010[ ). 

Real-space measurements and associated structure factors. We concluded above that 
a growing peak in Xi{t) 'directly' reveals the growth of a dynamic correlation length- 
scale as the glass transition is approached. This can only be correct if the assump- 
tions made above for the scaling form of G4(r;t) are correct. Dynamic lengthscales 
should in principle be obtained by direct measurements of a spatial correlation func- 
tion ( Doliwa and Heuer, 2000||Lacevic et al, 2002 1. 



However, in contrast to Xi(t), detailed measurements of Gi(r;t) are technically 
more challenging as dynamic correlations must now be resolved in space over large 
distances with a very high precision, and so there is much less data to draw on. From 
the point of view of numerical simulations where many measurements of Xi were 
reported, the main limitation to properly measure £4 is the system size. This might 
seem surprising as typical numbers extracted for the correlation length scale £4 are 
modest, growing, say, from 1 to at most 5-10 in most reports. However, this modest 
number hides the fact that the correlation function only decays to zero for distances 
r that are several times larger than £4. Given that G±{r;t) is accurately measured up 
to r = L/2 in a periodic system of linear size L, going to few £4 (say, five times), 
when £4^5 requires systems containing at least JV ~ L 3 ~ (2 x 5 x 5) 3 ~ 10 5 
particles in three dimensions, assuming the density is near p « 1. Such large system 
sizes are not easily studied at low temperatures when relaxation times get very large, 
even with present day computers. However, these studies arc of vital importance in 
that they allow the dynamical length scale £4(4) to be measured directly. Moreover, 
insights from such studies can then be used when inferring the behaviour of from 
measurements of Xi{~t)- 

Some representative data are shown in Fig. 10.71 They are obtained for a lattice gas 
model with kinetic constraints, where measurements of Gi(r;t) are somewhat easier 
than in molecular dynamics simulations. As discussed above, the idea is that for large 
r, G A {r;t) ~ A{t)r^F{r/^{t)) which then yields 

X 4(t) * A(t)Ut) d ~ P - (0-14) 

Often, one estimates the prcfactor A(t) to be equal to Gi{Q,t), which is simply the 
variance of the local quantity c(r; 0, t). However, the accuracy of this estimate is hard 
to assess without explicit evaluation of G±{r, t). In the example shown in Fig. 10.71 
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Fig. 0.7 Four-point correlation function Gi(r;t = r a ) measured in computer simulations 
of the Kob- Andersen kinetically constrained lattice gas in three dimensions ( |Berthie~2 003 ) . 
The dynamics slows down when density p increases, and the slower spatial decay of G4 directly 
reveals increasingly longer ranged dynamic correlations accompanying the glass transition. 

for instance the scaling between % 4 an d £4 in Eq. (|0.14[) is well obeyed, and careful 
examination of G 4 (r;i) suggests that p ra 1 and A is indeed a constant of order 1. 

While this is a subtle situation which requires each case to be considered indi- 
vidually, the work in this domain is broadly consistent with Xi(t)/Gi{Q,t) represent- 
ing the number of particles involved in heterogeneous relaxation. Note that these 
issues will also be relevant for discussion of other correlations and susceptibilities 
in later sections. Therefore, truly 'direct' measurements indeed confirm that the in- 
crease of the peak of x 4 (t) corresponds, as expected, to a growing dynamic lengthscale 
£ 4 (i) dDoliwa and Heuer, 2000|Lacevic et al, 2002|Bennemann et al, 1999|Berthier, 2004| 
|Berthier et al, 2007a[ ). As a result of subtleties related to the difference between four- 
point correlations in spin glasses and structural glasses, an early study of four-point 
functions ( |Dasgupta et al, 1991| drew an opposite conclusion, but the data of that 
study are in fact consistent with the now-established picture of a growing length scale. 

We also note, in passing, that the power p may have more than one interpretation. 
Typically, one assumes that a typical cluster has size £4^) and contains £ 4 (£) d_p parti- 
cles, so that d — p is interpreted as a fractal dimension. However, an alternative would 
be that clusters are all compact, but that the distribution of their sizes is rather broad. 
This uncertainty reflects the fact that four-point functions involve averages over many 
clusters, so that they do not resolve details of cluster structure. 

Instead of direct inspection of G 4 (r; t), it is often convenient to consider its Fourier 
transform S4(k;t) = J d d re lk r G 4 (r; t) ( |Bennemann et al, 1999[|Donati et al, 1999 



Glotzer et al, 2000||Lacevic et al, 2"002"l|Berthier, 2004[|Yamamoto and Onuki, 1998a 
Yamamoto and Onuki, 1998fe|Flenner and Szamel, 2010 ). In particular, this allows data 



for different times and different temperatures to be combined into a scaling anal- 
ysis that can yield the temperature dependence of £4, leaving an uncertain prefac- 
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tor. This approach was taken for example in Refs. ( |Yamamoto and Onuki, 1998 &[ 
|Berthier et ai, 2007a[|Flenner and Szamel, 2010[ ). Typically, simulation data result in 
length scales between 1 and 5 diameters. 

So far, we have considered circularly-averaged G4<V; t) and S±(k\ t). However, if the 
probe function Ci(t, 0) is anisotropic, as in Eq. (|0. 13[) . four-point functions G^r; t) then 
also depend on the angle between the separation vector r and the probe vector q. Sev- 
eral papers ( [Weeks et ai, 2007|Flenner and Szamel, 2009|Doliwa and Heuer, 2000[ ) have 
investigated this issue and found that motion in longitudinal directions is indeed more 
strongly correlated than motion in transverse directions, such that G^{r;t) is truly 
an anisotropic function. These findings add further to the difficulty of extracting the 
length scale £4 from direct measurements of four-point structure factors. 

0.2.5 Experimental estimates of multi-point susceptibilities 

Although readily accessible in numerical simulations, the fluctuations of C(t,0) that 
give access to X4 W are m general very small and impossible to measure directly in ex- 
periments, except when the range of the dynamic correlation is macroscopic, as in gran- 
ular materials ( |Marty and Dauchot, 2005] ) or in soft glassy materials where it can reach 
the micrometer and even millimetre range ( [Mayer et ai, 2004|Duri and Cipelletti, 2006[ ). 
To access Xi(t) m molecular liquids, one should perform time- resolved dynamic mea- 
surements probing very small volumes, with a linear size of the order of a few nanome- 
ters. Although possible, such experiments remain to be performed with the required 
accuracy. 

It was recently realized that simpler alternative procedures exist. The central idea 
underpinning these results is that induced dynamic fluctuations are in general more 
easily accessible than spontaneous ones, and both types of fluctuations can be related 
to one another by fluctuation-dissipation theorems. The physical motivation is that 
while four-point correlations offer a direct probe of the dynamic heterogeneities, other 
multi-point correlation functions give very useful information about the microscopic 
mechanisms leading to these heterogeneities. 

For example, one might expect that the slow part of a local enthalpy (or energy, 
density) fluctuation Sh x (t = 0) at position x and time t = triggers or eases the dy- 
namics in its surroundings, leading to a systematic correlation between Sh x (t = 0) and 
c(x + r; i, 0). This physical intuition suggests the definition of a family of three-point 
correlation functions that relate thermodynamic or structural fluctuations to dynam- 
ical ones. Interestingly, and crucially, some of these three-point correlations arc both 
experimentally accessible and give bounds or approximations to the four-point dy- 
namic correlations( |Berthier et ai, 20 05, Bcrthier et ai, 2007a, Bcrt hier et ai, 200"76| . 

Based on this insight, one may obtain a lower bound on Xi{t) using the Cauchy- 
Schwartz inequality (6H(0)5C(t, 0)) 2 < (SH(0) 2 ) (SC(t,0) 2 ), where H(t) denotes the 
enthalpy at time t, and SX = X — (X) denotes the fluctuating part of the observable 
X. By using a fluctuation-dissipation relation valid when the energy is conserved by 
the dynamics, the previous inequality can be rewritten as ( |Berthier et a/., 2005[ ): 

X4(t)>—{XT(t)}\ (0.15) 
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where the multi-point response function Xr(i) is defined by 

d(C(t,0)) 



XT® 



dT 



(6H(0)8C(t,0)}, (0.16) 



N,P 



and cp is the specific heat per particle (at constant pressure). 

In this way, the experimentally accessible response xt(^) which quantifies the sensi- 
tivity of average correlation functions (C(t, 0)) to an infinitesimal temperature change, 
can be used in Eq. (|0.15[) to yield a lower bound on Xi{i)- From Eq. (|0.16[) . it is clear 
that xt is directly related to the covariance of enthalpy and dynamic fluctuations, 
and thus captures the part of dynamic heterogeneity which is triggered by enthalpy 
fluctuations. 



Detailed numerical simulations and theoretical arguments (Berthier et al., 2007a 



Berthier et al., 20076 1 strongly suggest that the right hand side of (|0. 15p actually pro- 
vides a good estimate of X^W) m supercooled liquids, and not just a lower bound. 
Similar estimates exist considering density as a perturbing field instead of the tem- 
perature. These are useful when considering colloidal or granular materials where 
the glass transition is mostly controlled by the packing fraction. The quality of the 
corresponding lower bound was tested experimentally on granular packings close to 
the jamming transition (Lcchcnault et al., 2008 6| , and numerically for colloidal hard 



spheres ( |Brambilla et a/. ,12 009) 



Using this method, Dalle-Ferrier et al. ( |Dalle-Ferrier et al., 2007] ) have been able 
to estimate the evolution of the peak value of Xi f° r many different glass-formers in 
the entire supercooled regime. In Fig. 10.81 we show some of these results as a func- 
tion of the relaxation timescale. The value on the y-axis, a bound on the peak of 
is a proxy for the number of molecules, N CO it,4, in a cluster of mobile or immobile 
particles. As discussed briefly above, Xi(t) is expected to be equal to iVcorr,4, up to 
a proportionality constant A(t) which is not known from experiments, probably ex- 
plaining why the high temperature values of iV corr) 4 are smaller than one. Figure 10.81 
also indicates that N COTI< 4 grows faster when r Q is not very large, close to the onset 
of slow dynamics, and a power law relationship between -/V corr 4 and r Q is good in this 
regime (T a /r < 10 4 ). The growth of iV corri 4 with r a becomes much slower closer to the 
glass transition temperature T g , where a change of 6 decades in time corresponds to a 
mere increase of a factor about 4 of N corr ^, suggesting logarithmic rather than power 
law growth of dynamic correlations. A similar crossover towards a very slow growth 
of dynamic correlations is reported for colloidal hard spheres ( Brambil la et al., 2009[ ) 
and model glasses ( |Berthier et al., 2007a[ ), and is observed in numerical simulations 
even if the dynamic lengthscale £4 is directly estimated ( |Flenner and Szamel, 2010] ) . 
The consequences of such an effect for theories of the glass transition are discussed 
below. Bearing in mind all the caveats discussed above (unknown prefactors, quality 
of the bound, etc.), the experimental data compiled in Fig. 10.81 do appear to confirm 
that dynamic fluctuations and correlation lengthscales do grow when the molecular 
liquids approach their glass transitions. 
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Fig. 0.8 Dynamic scaling relation between the number of dynamically correlated particles, 
iVcorr,4, and relaxation timescale, T a , for a number of glass- formers fDallc-Fcrri er et al, 2007] ), 
determined using the bound provided by Eq. (|0.15[) . For all systems, dynamic fluc- 
tuations increase when the glass transition is approached. The full line through the 
data ( Dalle- Fc rrier et al., 2007| suggests a crossover from algebraic, N CO ti,a ~ t„, to loga- 
rithmic, A^corr.4 ~ exp(T„), growth of dynamic correlations with increasing r a . 

0.2.6 Four-point susceptibilities: some caveats 

The above story about the four-point susceptibility looks quite enticing. We have 
essentially argued that dynamical heterogeneity should be quantified by the spatial 
correlations of the mobility. This correlation function is a priori hard to measure in 
molecular glasses, but a divine surprise occurs: using rather trivial mathematics, its 
spatial integral is found to be bounded from below by a quantity that is much easier 
to measure. Is this too good to be true? What is the physics underpinning this 'easy' 
bound? 

The answer is that four-point correlation functions pick up a contribution that 
depends both on the statistical ensemble used (i.e. NVE vs. NVT) and on the dy- 
namics (i.e. Brownian vs. Newtonian). Using general scaling arguments based on a 
dominant length scale £4, one can show that in Fourier space the four-point correla- 
tion function has the following structure (in the k — > limit) (Bcrthie r et al., 20074 
|Berthier et al., 2007 

S 4 (k;r a ) = £fiTi(fc£ 4 ) + 9(k) kl^iM*)] 2 , (0.17) 



where s is a certain exponent (related to p above through 2s = d — p) and H\p, are 
certain scaling functions which behave similarly at small and large arguments. The 



Theoretical discussion 19 



wavevector dependent function g(k) carries the dependence on the statistical ensem- 
ble and dynamics. In particular, g{k = 0) = when all conserved quantities are 
strictly fixed, as in the NVE ensemble^. Similarly, g(k) is very different for Brown- 
ian dynamics (for which energy is not conserved) or for Newtonian dynamics (that 
conserves energy). So one has to face the uneasy truth that X4,(t) = Si(k = 0;t) de- 
pends on many microscopic and macroscopic details, although physically these should 
not affect the dynamic correlation length £4. It is definitely not so easy to directly 
relate Xi t° a dynamic correlation length. Intuitively, conserved variables play a role 
in "transmitting" the information about mobility from one region to another. 

From the above expression, one sees that ^(fc;^) mixes up two contributions that 
one would like to disentangle so as to extract the relevant scaling contribution from the 
first term only. It turns out that this first term is proportional to a three-point response 
function, that measures the change of the dynamics induced by a perturbation some 
distance away, and that we will discuss in more detail in section 10.41 It is this three- 
point response, not G4, that is the fundamental object carrying information about 
dynamical heterogeneities, and from which G4 is constructed. The space integral of 
the three-point response function is £|i/i(0); physically it represents the response of 
the dynamics to a uniform shift of an external parameter, such as the temperature or 
the density. Hence Xt(£) as defined in the previous section is a three-point function at 
zero wavevector. This explains the physical nature of the lower bound on Xi{t)'- due to 
the contribution of the energy as a conserved quantity, Xi has two contributions: one 
proportional to xt and one proportional to Xt> the latter being precisely the lower 
bound of the previous section. 

This discussion leads to the following caveats, that we alluded to above: a) the 
identification of a correlation volume from Xi(t) alone is not warranted in general. 
The information contained in S4 (fc; r Q ) is needed to unambiguously relate the growth 
of Xi(t) t° a growing lengthscalep b) extracting information from S4 can be difficult 
due to interference effects between the two terms and c) three-point response functions 
are the fundamental building bricks for dynamic correlations, and are not soiled with 
problems related to conserved variables or statistical ensembles. 

0.3 Theoretical discussion 

0.3.1 Recent progress based on four-point functions 

In the previous section, we have summarised some of the properties of four-point 
functions, their advantages for calculating the extent of dynamical heterogeneity, and 
some direct and indirect measurements of these quantities. 

There are many subtleties associated with these measurements, but the same broad 
picture is observed in a variety of systems and is robust to the precise measurement 
used. Essentially, as relaxation times increase, the four-point susceptibility increases, 

2 Actually, fast degrees of freedom can give subleading contributions in some cases. 

3 Related to this point, it is worth mentioning the case of purely Arrhenius systems, which are 
considered to be non-cooperative systems. Still, the lower bound based on xr proves that Xi( T a ) 
diverges in these systems as least as T~ 2 as the temperature goes to zero ( Dalle- Fcr rier et ah, 2007] l. 
The physical interpretation of such an apparent growth of the range of dynamical correlations is still 
unclear. 
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suggesting the presence of a growing dynamic correlation length. Where the real-space 
function Gn{r\t) can also be measured, this confirms more directly the increase of 
such a length scale, typically in the range between 1 and 10 molecular diameters. 
The fundamental, unavoidable conclusion seems to be that glassy behaviour is not a 
purely local 'caging' of particles by their neighbours, but indeed a genuine collective 
phenomenon. 

Having established the existence of a growing correlation length, several questions 
arise. From a theoretical perspective, we are familiar with the idea, borrowed from 
equilibrium critical phenomena, that when correlation length scales get large, micro- 
scopic features of the system become unimportant, and 'universal' behaviour emerges. 
Whether realistic glassy systems have length scales that are large enough for such a 
universal description remains unclear. Many analyses in this spirit have nevertheless 
been attempted, as we shall discuss shortly. It is likely that to in order to reach a good 
quantitative agreement a careful treatment of pre-asymptotic effects will have to be 
performed. 

A second fundamental point concerns the microscopic mechanisms that give rise 
to the correlations revealed by four-point functions. Many model systems can demon- 
strate the presence of increasing time scales, coupled with increasing susceptibilities 
Xi{t) and length scales £4. Predictions from different theoretical frameworks of the 
form of four-point functions are discussed in the next section, and we evaluate some 
of the theories in the light of existing results. 

0.3.2 Models of the glass transition and their predictions of dynamic 
heterogeneity 

We now turn to perhaps the most fundamental question in this area: what are the 
dominant mechanisms by which structural relaxation takes place in glassy materials? 
We give a quick survey of the dominant pictures for molecular glasses, their predictions 
for four-point functions, and the extent to which these are borne out. 

Mode-coupling theory. The mode-coupling theory (MCT) of the glass transition flGotze and Sjogren, 1 995 ) 
was historically derived from liquid state theory. Starting from exact microscopic equa- 
tions of motion for the density field in a liquid, several uncontrolled approximations are 
then performed to yield a closed set of dynamical equations for intermediate scattering 
functions. These equations give rise to a dynamic singularity at some finite temper- 
ature, T c , where relaxation times diverge in an algebraic manner. Additionally, very 
precise quantitative predictions can be made about the specific form of intermediate 
scattering functions, suggesting a very rich behaviour of time correlation functions, 
which do resemble the behaviour observed experimentally and reported in Fig. 10.21 As 
is well-known these predictions only apply over a modest time window of about 2-3 
decades in the moderately supercooled regime, but dramatically break down nearer to 
the glass transition QGotze, 2008[ ). 

Another dramatic failure of the traditional formulation of the theory, more rel- 
evant to the present contribution, is its inability to accurately predict the shape of 
the van-Hove distribution function described above, the resulting wavevector depen- 
dence of the self- intermediate scattering function (especially at low wavector) , and the 
corresponding decoupling between self-diffusion constant and the viscosity. 



Theoretical discussion 21 



Following this historical route, therefore, it is not obvious whether the MCT dy- 
namic singularity is accompanied by non-trivial dynamic fluctuations. This also means 
that the theory is not easily interpreted in physical terms. Recently, mode-coupling the- 
ory was reformulated in such a way that both these issues were greatly clarified. Using a 
field-theoretic formulation, it is possible to perform consistent mode-coupling approx- 
imations to get analytical predictions for both averaged two-time dynamic correlation 
functions and for the dynamic fluctuations around the averaged behaviour, i.e. for 
Xi{t) ( [Franz and Parisi, 2000l|Berthier et ai, 20076||Biroli and Bouchaud, 20040 . In a 



subsequent recent move, 'inhomogeneous mode-coupling' predictions for the shape and 
scaling form of four-point spatial correlation functions Gi(r;t) and its Fourier trans- 
form Si(k;t) were finally obtained ( |Biroli et al, 20060 . Thus, overall, mode-coupling 
theory is now able to make an impressive set of very detailed predictions for a very 
large family of spatio-temporal correlation functions for any given liquid, starting from 
the form of the microscopic interaction between the particles. 

A few numerical simulations have been presented to test these new predictions. 
First, the temporal evolution of the four-point susceptibility Xi(t) was compared to 
mode-coupling predictions. Just as time-correlation functions decay within MCT in a 
two-step process similar to the data presented in Fig. 10.21 Xi{t) is predicted to grow 
with time with two distinct power laws, \4 ~ t a an d Xi ~ t b m the time regimes respec- 
tively corresponding to the approach to, and departure from, the plateau. These two 
power law regimes have been successfully identified in numerical work, with numerical 
values for the exponents a and b that are in 'reasonable' agreement with numbers 
predicted by MCT ( |Berthier et ai, 20076[|Berthier, 2007[|Berthier and Kob, 20070 . 

The peak of the four-point susceptibility is predicted to diverge algebraically at the 
critical temperature. This prediction was observed numerically to hold over a similar 
(restricted) temperature window as for the averaged relaxation time r Q itself. This 
finding implies that the peak of Xii when plotted as a function of r Q , follows a power 
law scaling, Xi ~ t qi where z is predicted to be a non-universal critical exponent. 
Returning to the data compilation in Fig. 10.81 we remark that the data obtained in 
the moderately supercooled regime do indeed approximately follow an initial growth 
which is consistent with the MCT prediction, while clearly breaking down at lower 
temperatures. 

Finally, predictions for the detailed shape of four-point correlation functions, in 
particular for Si(k; t), were recently confronted to numerical results, with inconclusive 



results. While a first paper (Stein and Andersen, 20081 reports excellent agreement 



with MCT predictions both for the wavevector dependence of Si(k;t) and its evo- 
lution with temperature, a more recent report (Karma kar et at, 20090 claims that 
disagreements with theoretical predictions arise when larger system sizes are included 
in the numerical analysis. This ongoing debate illustrates the fact mentioned above 
that even for the modest correlation length scales characterising relaxation in super- 
cooled liquids, very large system sizes are needed to unambiguously and accurately 
measure four-point functions. Clearly, more work is needed to clarify the status of the 
large body of MCT predictions regarding four-point functions. 



Facilitation picture and kinctically constrained models. In the facilitation picture of 
supercooled liquids, structural relaxation is thought to originate from propagation of 
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localised excitations, or 'defects', throughout the system flGlarum, 1960D - The idea 
is that the local structure at position x changes rapidly when a defect visits the 
neighbourhood of x. Then, one makes the further hypothesis that defects are sparse, 
and most of the system is in fact immobile, such that defects can be considered as 
uncorrelated, independent objects. In practice, very little is known about the nature, 
origin, or even the existence of such defects in real liquids. 

Nevertheless, it is clear within this picture that dynamics is highly heterogeneous 
in space, and temporally intermittent. Physically, one expects a dynamic correlation 
lengthscale £4 to emerge, which should basically correspond to the linear size of the 
region explored independently by a given defect. This also means that the time depen- 
dence of £4(4), and thus of Xi{t) can directly be connected, in this view, to how fast 
the defects move. In the simplest approximation where defects are diffusive objects, 
one would expect ^(t) ~ \ft for t < r Q . 

In recent years, these ideas have been pursued quite extensively, based on the pro- 



posal (Garrahan and Chandler, 2003) that facilitation is essential for explaining the 
heterogeneous dynamics of supercooled liquids. The theory has been developed primar- 
ily through studies of systems called kinetically constrained models fRitort and Sollich, 2003[ ). 
Numerous distinct lattice models belong to this family, which can be distinguished by 
the set of microscopic rules governing the dynamics of the localised defects, the exis- 
tence or absence of conservation laws, the topology of the lattice, etc. The simplest 
models, such as the (one-spin facilitated) Fredrickson-Andersen (FA) model, in fact 
reproduce nearly exactly the scaling relations mentioned above for four-point func- 
tions ( |Toninelli et al, 2005|Chandler et al, 2006|Whitelam et al, 2004|Whitelam et al, 2005[ ). 

Given the large number of distinct models, and the fact that models are postulated 
instead of being derived as approximate, coarse-grained representations of liquids, it is 
not clear how one should compare their behaviour to numerical results obtained for re- 



alistic liquid models. This issue is discussed in a recent review ( Chandler and Garrahan, 2010 ) 
Thus, it is perhaps better to interpret this diversity as being suggestive of the differ- 
ent types of behaviour one can possibly encounter in liquids. On the other hand, from 
the theoretical point of view, having well-defined, relatively simple, statistical models 
defined on the lattice is very appealing as very many detailed and quantitative re- 
sults can be obtained by exploiting tools from statistical mechanics. Thus, kinetically 
constrained models can also be viewed as 'toy supercooled liquids'. 

In this regard, the study of the dynamically heterogeneous behaviour of kinet- 
ically constrained models has been a very active field of research in recent years. 
Many models have been investigated and a large number of time correlation functions 
(two-point, four-point, persistence functions) have been analysed, suggesting possible 



behaviours for dynamic susceptibilities (Berthier et al, 2007a Chandler et al, 20061 



or decoupling phenomena flJung et al, 2004 Pan et al., 2005 1 . We refer to the chapter 



by Sollich, Toninclli, and Garrahan for further details and references on this topic. 

At the qualitative level, it is obvious that all models are characterised by rapidly 
growing time scales and length scales, and are thus interesting models to study dy- 
namic heterogeneity. However, models with diffusive point defects (like the simplest 
of Fredrickson-Andersen models), do not compare well with the real liquids that have 
been studied so far. In three dimensions, they predict simple exponential relaxation 
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and no decoupling phenomena ( |Jack et al, 2006] ) . The dependence of X4 on time and 
temperature are also characterised by scaling laws that have not been observed in 
numerical and experimental results ( jTonineili et al, STO OSiBcr thier et al, 2005] ). How- 
ever, the Arrhenius scaling of their relaxation time indicates a relation between these 
models and 'strong' liquids (Garrah an and Chandler, 2 003, Whitcla m et al, 2005 ) and 



there are comparatively few results on dynamical heterogeneity for such materials. In 
particular, while such models predict rather large dynamical length scales in strong ma- 
terials QGarrahan and C handler, 2003, Bcrth ier and Garrahan, 2005] ), these have not 
yet been observed. 

On the other hand, more complicated models where defects move sub-diffusively or 
cooperatively seem to be more appropriate representations of 'fragile' liquids, which 
have a non-Arrhenius scaling of relaxation time with temperature. Such kinetically 
constrained models exhibit stretched exponential relaxation as in Fig. 10.21 decou- 
pling phenomena QJung et al, 2004| similar to the results in Fig. 10.41 realistic form of 
four-point structure factors as in Fig. 10.71 or dynamic length scales which grow very 
slowly ( jGarrahan and C handler, 2003), in qualitative agreement with the experimen- 



tal results shown in Fig. 10.8 

Adam-Gibbs and the mosaic picture. The idea that relaxation events in glasses are 
collective and involve the simultaneous motion of several particles dates back at least to 
Adam and Gibbs, who provided an argument to relate the size of these "cooperatively 
rearranging regions" (CRR) to the configurational entropy of the supercooled liquid. A 
lot of the work on dynamical heterogeneities is in fact motivated by the Adam-Gibbs 
picture and attempts to determine the size of these CRR ( |Binder and Kob, 2005[ ). 

The Adam-Gibbs picture was later put on more solid ground in the context of the 
Random First Order Transition (RFOT) Theory of glasses ( |Kirkpatrick et al, 1989} 
Lubchenko and Wolynes, 20 06). RFOT suggests that supercooled liquids be can be 



thought of as a mosaic of "glass nodules" or "glassites" with a spatial extension £* (T) 
limited by the configurational entropy. Regions of size smaller than £*(T) are ideal 
glasses: they cannot relax, even on very long time scales, because the number of states 
towards which the system can escape is too small to compete with the energy that 
blocks the system in a given favourable configuration. Regions of size greater than 
£* are liquid in the sense that they explore with time an exponentially large number 
of unrelated configurations, and all correlation functions go to zero. The relaxation 
time of the whole liquid is therefore the relaxation time of glassites of size £* . This 
relaxation occurs through collective activated events that sweep a region of size £*, 
which are the CRR regions of the Adam-Gibbs theory. The crucial assumption of 
RFOT is that thermodynamics alone fixes the value of £* , whereas the relaxation time 
r a involves the height of the activation barriers on scale £* , which is assumed to grow 
as a power-law, £*^ , where tp is a certain exponent ( |Bouchaud and Biroli, 2004| . 

Note that when an activated event takes place within a glassite of size I*, the 
boundary conditions of the nearby region changes. There is a substantial probability 
that this triggers, or facilitates, an activated event there as well, possibly inducing an 
"avalanche" process that extends over the dynamic correlation length scale £ > £* . The 
dynamics on length scales less than £* is, within RFOT, inherently cooperative, but 
the relation between the dynamic correlation length £, defined for example through 
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three- or four-point point correlation function and the mosaic length t* is at this stage 
an important open problem (see ( |Dalle-Ferrier et ai, 2007[[Capaccioli et ai, 2008[ )). 
If £ is of the order of a few glassite lengths I* , then one expects that X&{ T a) should 
grow as t to some power. Assuming activated scaling, lnr Q ~ finally leads to 
X4 ~ (In t q ) z , instead of a power-law relation predicted by MCT or by non-cooperative 
KCMs. The crossover towards this logarithmic behaviour is not incompatible with 
the data ( |Xia and Wolynes, 2000[ ), see section ID. 2.51 above. In fact, the details of the 
crossover between the MCT region and the RFOT region are still very mysterious 



(see ( |Biroli and Bouchaud, 2009| ) for a recent discussion), but some claims have been 
made about the evolution of the shape of the dynamically correlated regions, that 
should morph from stringy, fractal objects in the MCT region to compact blobs at 
lower temperatures ( |Stevenson et at, 2006[ ). It would be interesting to devise some 
experimental protocol to test these predictions. 

0.4 Beyond four-point functions: other tools to detect dynamical 
correlations 

So far, we have discussed how four-point functions can be used to estimate dynamical 
length scales, and we have stated that these are typically found to be in the range 1-10 
molecular diameters. However, we have also noted that (a) the four-point susceptibility 
estimate of the dynamical correlation volume may lead to erroneous results (see section 
I0.2.6j) and (b) a variety of different theories of the glass transition are broadly consistent 
with the above estimate. 

To make further progress, it seems that more adapted and discriminating observ- 
ables will be required. In fact, there are a wealth of methods that have been used 
to characterise dynamical heterogeneity, of which we discuss just a few, and we refer 
to other chapters in this book for details. Here, we mainly emphasise the questions 
that can be addressed by different methods, and give an overview of the relationships 
between some of the methods that have been developed in different contexts. 

0.4.1 Non-linear susceptibilities 

In standard critical phenomena, diverging two-point correlations lead to singular lin- 
ear responses. It is therefore quite natural to conjecture that increasing dynamical 
correlations should also lead to anomalous responses of some kind. Spin glass theory 
provides, again, an interesting insight. As discussed in section 10.2.41 the spin glass 
transition (at zero external field) is signaled by the divergence of the four-point static 
correlation function xsg- It can be easily established that close to the transtion, xsg 
is related to the third-order non- linear magnetic response at zero frequency \3 ( w = 0) 
dFischer and Hertz, 1991] ): 

/ n x Xsg - 2/3 ,„ 1Q , 

X3(w = 0) = "W ( } 

Thus, although linear responses are blind to the development of spin glass long range 
order, the non-linear magnetic response is not. Actually, it diverges at the transition 
and, hence, is a direct experimental probe, contrary to xsg, which can instead only 
be measured in numerical simulations. 
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The analogy with spin glasses discussed in section 10.2.41 therefore suggests that 
glasses should also display increasingly non-linear responses approaching the glass 



transition, as first argued in (Bouchaud and Biroli, 2005). Theoretically, this can be 



substantiated by some general scaling arguments and by a mode-coupling calculation. 
These are described in ( |Tarzia et gl., 2010[ ); we will just briefly summarize them in the 
following. The starting point is to rewrite the generic third order non-linear response 
X3 (t) in terms of the second order change R 2 of the linear response R: 

X.W = ^^» ^)^^ ^0^W) = (0-19) 

/ dt 1 dt 2 dt 3 E(t 1 )-^^—E(t 2 )E(t 3 )=[ dt 1 E(t 1 )R 2 (t,t 1 ) 
J-oc dE(t 2 )dE(t 2 ) J.oo 

Note that we will focus on the dielectric non-linear response (so P is the electric 
polarization) but the generalization to other perturbing field is straightforward^. It is 
easy to understand, at least at low frequency, why R 2 and therefore the non-linear 
susceptibility have a singular behavior. In fact, within an adiabatic approximation, 
one finds that the linear response in the steady state created by a slowly alternating 
field is: 

R eq (t-t',Ecos(ujt)), (0.20) 

where R eq (t — t',E) is the equilibrium response function with a static field E. Since 
we are interested in the small E behavior, we can expand the above expression up to 
second order in E, this yields: 



R 2 (t,h) =R eg (T,Ecos(ut))-R (T) 



E 2 cos 2 '(ut) d 2 R eq {t-t 1 ,E) 



dE 2 



, (0.21) 

E=0 



where Ro(t) is the unperturbed equilibrium response function, and the derivative is 
computed with respect to a constant external field. The second term is expected to 
give a singular contribution because close to the glass transition a small applied field 
E is roughly equivalent to a shift of the order of E 2 of the glass transition temperature 
(see below), and T a significantly varies when the temperature is changed by a small 
amount close to T g . More precisely, by taking into account the E — s- —E symmetry 
one can rewrite d 2 R eq /dE 2 as 2dR eq /dO where 9 — E 2 . Using the time temperature 
superposition, one finds: 



dR eq (r, E) 



g« ' ■ (0-22) 

E=0 01 



where k = d J T a ^® . Using that R eq (r, 0) is the Fourier transform of the linear suscep- 
tibility Xi(w) an< i P nl gg m g th° previous expressions in (|0.19j) (and after some algebra 
detailed in ( |Tarzia et ai, 2010] )) one finds the following result: 

4 Depending on the perturbing field, the symmetry E — > —E will hold or not. In the latter case, 
the first non-linear response is the quadratic one, in the former the quadratic vanishes by symmetry 
and one has to focus on the third order one. 
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Fig. 0.9 Sketch of log |x3(w)| as a function of logcj, showing different frequency regimes and 
crossovers ( |Tarzia et al, 2010| : u)T a <C 1, uT a ~ 1, Tg jr a <C 0JTp <C 1 (e = T — T c ), o;t/3 2> 1, 
loto ~ 1. Note that the low frequency limit is non zero but much smaller than the peak value 
for T close to T c . 

X3[u) « « — , (0.23) 

which is expected to hold at low enough frequency, at least when the deviations from 
time temperature superposition are weak, k is expected to be a slowly varying func- 
tion of temperature, a constant in first approximation ( |Tarzia et al., 2010[ ). In this 
expression, 9%i(w)/9T is akin to the three-point susceptibility xt defined in section 
10.2.51 Thus, Eq. (|0.23[) is an important result since it establishes a relationship with 
the linear dynamical responses that have been used to evaluate dynamical correla- 
tions, and it also proves that supercooled liquids should respond in an increasingly 
non-linear way approaching the glass transition since, as we have discussed before, xt 
and therefore increase approaching the glass transition. 

The above general heuristic arguments can be supplemented by more microscopic 
ones based on MCT, which provides quantitative predictions on the critical behaviour 
of X3- Although the corresponding results are restricted to the small temperature 
regime where MCT is believed to apply, they are nevertheless guidelines for the general 
behaviour of X3- We sketch the evolution of the absolute value of \3 with frequency 
in Fig. MM 

• In the cv-regime, i.e. uj ~ 1/tc, - e^ 2a+1 / 2b /T , the absolute value of xsM grows 
with decreasing u and reaches its maximum, of height of order 1/e, after which it 
decreases as uj~ b Tp at large ui. In this regime, one has the scaling form: X3( w ) = 

• At the crossover between the early a-regime and late /3-regime the absolute value 
of X3( w ) is °f order l/^/e- 
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• In the /3-regime, i.e. 10 ~ ~ e 1 / 2 "/^, the absolute value of X3( w ) decreases 
as Lo~ b Tp at small lo and as oj~ a Tp at large w. In this regime one has a scaling 
form x 3 (w) = 

Exponents a and b are the well known critical exponents of MCT introduced above; 
To is a microscopic relaxation time and e = (T — T c )/T c is the distance from the mode- 
coupling critical temperature T c . We remark that the existence of the peak and the de- 
crease at low frequency is a non trivial prediction since it is in contrast with the (trivial) 
non- linear response of uncorrelated Brownian dipoles (Dejardi n and Kalmykov, 1999[ ) 
and with spin glasses. The decrease with an exponent three at high frequency sketched 
in Fig. IO.9l is instead trivial and present also for independent dipoles (Dojardi n and Kalmykov, 1999 1 



Recent experiments of third-order non-linear dielectric responses of supercooled 
glycerol have indeed shown, for the first time, that these theoretical expectations are 
qualitatively correct flCrauste-Thibierge et al., 2010[ ) . We refer to the chapter by Alba- 
Simionesco et al. for a presentation of these experimental results and their discussion. 

We conclude this section by emphasising that perturbing fields, other than electric 
fields, are expected to lead to similar results. Studying non-linear responses seems to 
us a very promising route to follow in order to probe the glassy state in a new way. 
A particularly interesting case worth studying corresponds to non-linear mechanical 
responses of colloidal glassy liquids. In this case, the values of the perturbing field that 
affect the sample are of the order of Pa, thus much smaller than the ones affecting 
the measuring apparatus, which arc of the order of GPa. Thus, within a very good 
approximation, the only non-linear output signal is from the sample itself. This is not 
the case for dielectric measurements, which are therefore very difficult since one has 
to be able to filter out the trivial non-linear part due to the amplifiers, etc. This was 



the main difficulty in the experiment reported in ( Crauste-Thibierge et al., 2010). 



0.4.2 Inhomogeneous dynamical susceptibilities 

We have seen in section 10.2. 51 that the variation of a dynamical correlator with respect 
to an external parameter (e.g., xt) is a way to obtain estimates of the number of dy- 
namically correlated particles. As a natural generalisation, one can study the variation 
of a local correlator, which measures the dynamics around the position x, induced by 
a perturbation at certain other point z. By summing over all z one obtains again a 
global dynamical response such as \t, since this then corresponds to computing the 
variation with respect to a uniform shift of the external parameter. 

This new, spatially dependent, dynamical response function is akin to G4 and allows 
one to probe the spatial structure of dynamic heterogeneity and to measure directly 
a dynamical correlation length £. The physical reason is that spontaneous dynamical 
fluctuations measured by the 4-point function and induced dynamical fluctuations 
measured by this new type of response function are intimately related. Accelerating 
or slowing down the dynamics at one given point (by adding an external potential) 
must perturb the dynamics over a length scale £ if the dynamics are indeed correlated 
over this distance. 

Let us define more precisely this new dynamical response. Consider the change 
in the local dynamical structure factor F(x, y, t) due to an extra, spatially vary- 
ing, external potential U (z) . Note that this observable can always be decomposed in 
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Fourier modes: 5J ^§f \ u=0 = J dkdq e - lq ( x - y ) +4k (y - z )xk(q^), where Xk(q,*) oc 

~ w'(k) k ' t '* l[/=o * s ^ ne response of the dynamical structure factor to a static exter- 
nal perturbation in Fourier space. For a perturbation localised at the origin, U (z) = 
UqS(z), one finds SF(q,y,i) = Uq J dke y Xk(q,i)- This susceptibility is also related 
to a 3-point density correlation function in the absence of the perturbation. It is very 
important to note the very different role played by q (the standard wavevector) and 
k (the wavevector of the perturbation): only the latter is sensitive to dynamic corre- 
lations n 

The susceptibility Xk(q, t) is interesting for experimental reasons, at least in col- 
loids, since it could be measured by using optical tweezer techniques. From a funda- 
mental point view, it provides a very useful way to characterise dynamic heterogene- 
ity since it is not affected by complications due to conservation laws and the type of 
dynamics, contrary to X4 an d G&. Therefore we expect that extracting spatial infor- 
mation and, especially, a precise estimate of £ should be cleaner by using Xk(q, t) (see 
the discussion in section l0.2.6j) . 

Finally, another advantage of Xk(q, t) is that precise quantitative predictions have 
been obtained within MCT by analytical arguments ( |Biroli et al., 2006] ), which were 
later confirmed by numerical analysis ( |Szamel and Flenner, 2009[ ) and complementary 



approaches (Szamel, 2008). The critical behaviour of Xk(q, t) approaching the MCT 
transition temperature is the following (we use the same MCT notation introduced 
previously): 

• In the /3-regime, i.e for times of the order of rg = e _1 / 2a , one finds 

* t(q ''' " TfTTS* ("S' i£l/ "°) ' (a24) 

where the proportionality constant depends on q. The scaling function gp{k 2 /y/e, tt 1 ^ 2a ) 
is regular for k = 0, thus implying that the k = value diverge as 1/y/e. For large 
values of u = ie 1 / 2 " one finds that gp{k 2 / y/e, te 1 / 2a ) equals E(Tk 2 / yfe)u b , with S 
a certain regular function. 

In the a-regime, i.e. for times of the order of t q = e - 1 / 2a - 1 / 2b ; onc finds 



Xk(q,*)= WP^ 9^ (^) 



yfe{yfe + Tk 2 )^ q \r L 

with S is the same function defined previously. It has the properties: S(0) 7^ 
and S(w 3> 1) ~ 1/v such that Xk behaves as fc~ 4 for large fee -1 / 4 , independently 
of e. Also, g a , q (u < 1) a u b , as to match the j3 regime, and g a (u 3> 1, k) — > 0. 

Note that the spatial scaling variable is k 2 jiT 1 ! 2 in both the a and the (3 regimes. 
The physical consequence is that there exists a unique diverging dynamic correlation 
length £ ~ v^lel" 1 / 4 that rules the response of the system to a space-dependent 
perturbation within MCT. The analysis of the early (3 regime where t <C |e| _1//2a 



5 We note that compared to Refs. flBcrthier et al., 2007a, Bi roli et al., 2006| l, the notations q and 
k have been inverted. 
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shows that this length in fact first increases as W 2 and then saturates at £. Interest- 
ingly, this suggests that while keeping a fixed extension £, the (fractal) geometrical 
structures carrying the dynamic correlations significantly "fatten" between t@ and r Q , 
where more compact structures are expected, as perhaps suggested by the results of 
dAppignanesi et ai, 2006[ ). 

Up to now, there are no simulations or experiments measuring a spatial dynamic 
response such as Xk(q, t). Hopefully, these will be performed in the future. As dis- 
cussed previously, we do believe that this new observable is a simpler and more direct 
measure of dynamical correlations than G4. Furthermore, quantitative results beyond 
scaling can be obtained within MCT. Thus one could consider comparing the MCT 
predictions for dynamical heterogeneities to numerical and experimental result in a 
stringent way, as it has been done for the intermediate scattering function, see e.g. 



(Kob and Andersen, 1995) 



0.4.3 Structure and dynamics: Is dynamic heterogeneity connected to 
the liquid structure? 

One of the most frequently asked questions in studies of dynamical heterogeneity 
is whether the observed fluctuations might be structural in origin. Such questions 
have attracted sustained interest. For example, in early numerical work on dynamic 
heterogeneity, immobile regions were discussed in connection with compositional fluc- 
tuations in fluid mixtures ( |Hurley and Harrowell, 1995 1. Thirteen years later, some 



form of local crystalline order is invoked to account for slow domains in numerical 
work (Kawasa ki et ai, 2007[ ). 

It should be noted that this chapter, and perhaps even this whole book about 
'dynamic heterogeneity' would not exist in this form if the question of the connec- 
tion between structure and dynamics had been satisfactorily answered. In that case, 
indeed, research would be dedicated to understanding the development of structural 
correlations at low temperatures in supercooled liquids, and to developing tools to 
measure, quantify and analyse such static features. 

A key advance in connecting structural properties to dynamical heterogeneity has 



been the development of the so-called 'isoconfigurational ensemble' QWidmer- Cooper et al., 2004 1 
In this approach, one calculates a traditional ensemble average in two stages. First, one 
averages the particles' velocities, keeping their initial positions fixed. (If the dynam- 
ics are stochastic, this step also contains an average over random noises.) Averaging 
a local dynamical observable such as c(r; t, 0) in this way, one arrives at an 'isocon- 
figurational average' (c(r; t, 0))j so , which still depends on the position r through the 
fixed initial particle positions. This average is therefore able to reveal the influence of 
the structure of the initial configuration on the dynamical behavior at that point. To 
return to a traditional ensemble average, one carries out an average over the initial 
particle positions in a second step. 

The right panel of Fig. 10.11 represents the spatial dependence of the isoconfig- 
urationally averaged single-particle mobility in a two-dimensional mixture of soft 
disks ( |Widmer-Cooper et ai, 2004[ ). The fact that this image is not uniform demon- 
strates that part of the dynamic heterogeneity has a structural origin. This raises two 
different questions. First, can one predict from structural measurements the pattern 
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produced by the isoconfigurational average in Fig. 10.11 ? Second, how much of the 'real' 
dynamic heterogeneity is actually preserved by the isoconfigurational average and has 
thus a genuine structural origin? 



Harrowcll and coworkers have provided detailed answers to the first question ( 


Widmer-Cooper et al., 2004 


|Widmer-Cooper and Harrowell, 2006|Widmer-Cooper and Harrowell, 2005|Wic 


mer-Coopcr and Harrowell, 2007D 



Statistical analysis of isoconfigurational ensembles has been very useful in assess- 
ing the statistical significance of correlations between mobility and the local energy, 
composition or free volume. They have recently made the point that strong corre- 
lations exist between vibrational properties of the liquid and isoconfigurational mo- 
bilities ( Widmer-Cooper et al., 2008[ ). They have also made vivid the distinction be- 
tween the existence of a statistical correlation between structural and dynamical fluc- 
tuations, and the much more demanding notion of a causal link between the two, 
that is, of a correlation that is strong enough that prediction of the mobility can be 



made based on a given structural information (Widmer-Cooper and Harrowell, 2006 
|Widmer-Cooper and H arrowel h~2005[ ). These two notions are very often confused in 
the dynamic heterogeneity literature. 

The example of enthalpy fluctuations is useful in this respect. The fact that the 
four-point susceptibility Xi{t) can be quantitatively estimated with good accuracy 
from a three-point susceptibility such as xt(*) °c (SH(0)SC(t,0)} provides evidence of 
a strong correlation between enthalpy and dynamic fluctuations. However, enthalpy 
fluctuations are not good predictors for dynamic heterogeneity, presumably because 
they contain short-ranged and short-lived fluctuations that do not correlate well with 
slow dynamics. Indeed, suitably filtered enthalpy fluctuations correlate very strongly 
with dynamic heterogeneity ( Matharoo et al., 2006[ ) 



We finally return to the second question: is dynamic heterogeneity truly captured 
by isoconfigurational averages, and thus does it fully originate from the structure? The 
response is more subtle than expected as it depends on which observable, and more 
precisely on which lengthscale, it is analysed. We mentioned above that dynamic het- 
erogeneity primarily revealed itself through the intermittent single particle dynamics 
(Fig. 10 . 3[) leading to broad distributions of single particle displacements with broad 
tails. These features almost completely disappear after the isoconfigurational average 
is performed (Berthier and Jack, 2007[ ). In other words, the distinction between mo- 
bile and immobile particles is mostly dynamical in nature, suggesting that the quest 
for a connection between the static and dynamic properties of glass-formers at the 
particle level is in vain. 

Nevertheless, mobility fluctuations do display interesting spatial correlations, as 
illustrated in Fig. 10.11 This suggests that the distinction between fast and slow do- 
mains remains consistent in the isoconfigurational ensemble. This observation can be 
quantified by measuring a 'restricted' four-point function 

XT (t) = iv((C(t, 0) 2 ) iso - (C(t, 0))? o ) • (0.26) 

\ / initial cond. 

While Xi(t) measures the total strength of dynamic heterogeneity, X4°(*) makes use of 
the isoconfigurational ensemble and first records the strength of dynamic heterogeneity 
at fixed initial conditions, the average over initial conditions being performed after- 
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wards. In the case where isoconfigurational mobility (and thus the image in Fig. 10. lj) 
is uniform, one has X4°W = X4Wi since the average over initial conditions is triv- 
ial in this case. More generally, a large contribution of X4°(*) to Xi{t) indicates 
that the dynamic fluctuations captured by Xi{t) are inherently dynamical in origin 
and do not originate in the liquid structure. Numerical measurements in molecular 
dynamics simulations indicate that the opposite is true and X4°M contributes less 
and less to x*(t) as temperature decreases ( |Berthier and Jack, 2007[ ). This suggests 
that the search of a causal link between structure and mobility does make sense, 
at least on large length scales. Interestingly, the vibrational properties investigated in 
Ref. ( Widmer-Cooper et al., 2008 ) as relevant structural indicators of dynamic hetero- 
geneity are a suitable candidate, since the vibrational spectrum is a collective property. 



0.4.4 Point to set correlations: Emergence of amorphous long range 
order? 

As discussed in the previous section, it is quite natural to ask what structural features 
(if any) might be responsible for the growth of dynamical correlations. One possibility 
is that actually there exists a static growing length that drives the increase of dy- 
namical correlations. As we discussed in the introduction, simple static correlations 
are rather featureless when approaching the glass transition. However a new length 



called the "point to set" length was recently introduced (Bouchaud and Biroli, 2004 



|Mezard and Montanari, 2006 1 . It is naturally devised to probe the growth of static 



amorphous long range order ( Bouchaud and Biroli, 2004 Mezard and Monta nari, 2006[ ) 



and has been shown to grow close to the glass transition ( Biroli et al, 2008 ). This is 
reviewed in the chapter by Semerjian and Franz. 

The basic idea is to measure how much boundary conditions affect the behaviour 
of the system, far away from the boundaries themselves. This is the usual way to test 
for the emergence of long range order in statistical mechanics. However, for standard 
phase transitions, the appropriate boundary conditions are known from the outset. 
For example, in the case of ferromagnetic transitions, one can fix the boundary spins 
mostly in the up direction and check whether this leads to a positive magnetisation for 
spins in the bulk. The difficulty in the case of glasses, for which one would like to test 
the presence of long range amorphous order, is that the boundary conditions one has to 
use look just as random as the amorphous configuration one wants to select. The way 
out is to let the system itself choose the boundary conditions: the procedure is to take 
an equilibrated configuration a, freeze all particles outside a cavity of radius 1Z and 
then recompute the thermodynamics for the particles inside the cavity, that now are 
subjected to a typical equilibrium boundary condition. One can then study a suitably 
defined average overlap q(1Z) between the new thermalized configurations at the centre 
of the cavity and the reference state a, as a function of TZ. The quantity q(TZ) is called 
a "point-to-set" correlation ( jMezard and Montanari, 2006[ ). The characteristic length 
scale over which q(TZ) drops to zero is called the point to set length. The increase of 
this length is a clear signal that the system is developing long range static order, and 
in the case of glasses, amorphous long range order. This point-to-set length is precisely 
the size of the glassitcs within the RFOT theory of glasses, see section lQ.3.21 
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This topic is discussed in detail by Franz and Semerjian, to which we refer for a 
presentation of the general theoretical and numerical results. Here we simply mention 
that this point to set length has been shown to grow in numerical simulations of su- 
percooled liquids (Biroli et ai, 20081. Furthermore, it was proved that it must diverge 
whenever the relaxation time does so ([Montanari and Semerjian, 2006[ ). 

Therefore, an important open question is whether the correlation length picked up 
by the dynamical correlators discussed above are actually just a consequence of hidden 
static correlations or if they are instead quite unrelated to them. In the first case, the 
study of dynamical correlation will still be very valuable because it provides an easier 
way to probe static correlations. Whatever is the correct answer, it would lead to a 
substantial progress in our understanding of the glass transition and would help us in 
pruning down the correct theory. We do hope that numerical and experimental studies 
will be devoted to this important problem in the future. 



0.4.5 Large deviations and space-time thermodynamics 

It is clear from its definition in Eq. (|0.7[) that the correlation G^r, t) is the covariance of 
the mobility c(r;0,i) at two nearby spatial points. Similarly, the three-point functions 
of the previous section are covariances of c(r;0,t) with the local energy, enthalpy or 
free volume. Of course, not all information about the mobility is contained in such 
covariances: one might consider higher moments of these functions or indeed the joint 
distribution of mobilities at all points. However, the inherent difficulty of characterising 
the distribution of an entire mobility field requires physical intuition in choosing which 
observable to measure. 

A natural first choice for such a scheme is to consider the fluctuations of the spa- 
tially averaged correlation function C (i, 0) , beyond the Gaussian level. Such measure- 
ments are possible in experiments such as those of Duri et al. ( Duri et ai, 2005[ ) (see 



Fig. IQ.lOj) and in computer simulations of a variety of models ( Chamon et ai, 2004[ ). 
They have typically been considered in out-of-equilibrium situations but this is not 
essential (see also below). Typically, such distributions are skewed and non-Gaussian, 
and it is natural to connect the asymmetry of the distribution with dynamical hetero- 
geneity. For example, even on time scales t much greater than the structural relaxation 
time, there is a substantial probability that regions of the system have persisted in 
an immobile state for all times between and t. This enhances the probability of 
observing a larger than average value for C(£,0). The opposite behaviour may occur 
on short times: while typical regions have not relaxed, co-operative motion in some 
regions enhances the probability of observing a smaller than average of C(t,0). 

In making this connection, it seems reasonable that regions where C(i,0) is large 
possess rather stable structure at the molecular level, while regions where C(t,0) is 
small correspond to relatively unstable local structure. This fact has recently been 
exploited in computational studies that probe trajectories where relaxation is much 
slower than average. To identify such trajectories, it is useful to define a measure of 
dynamical activity, for systems of N particles evolving over an observation time t bs- 
For example, one may take ( |Hedges et ai, 2 009) 
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Fig. 0.10 Distributions of dynamical observables. (Left) Distribution of the correlation func- 
tion in a coarsening foam ( |Duri et al, 2005| . The skewed distribution arises from rare trajec- 
tories with more motion than average. (Centre) The distribution of the 'dynamical action' e 
in a kinetically constrained model ( Mcrollc et al., 2005[ ) . The action e measures the amount 
of motion in a trajectory in the same spirit as the activity K defined in Eq. (I0,27[) . The 
distribution is left-skewed and non-convex, indicating a population of trajectories with low 
mobility. (Right) Distribution of the activity K in a model of Lennard- Jones particles, with 
a biasing field s in place, as discussed in the text. Two peaks are evident when the t b B is 
large, revealing the presence of two distinct dynamical phases ( |Hedges et al, 2009[ ). 
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where Vi(t) is the position of particle i at time t and the tj = jAt are equally spaced 
times. For large N and t bs, the distribution of K becomes sharply peaked about its 
average, (K). 

In general, for large N and i bs, one expects a the equilibrium distribution of K 
to have the form 

P{K) ~ exp[-Nt obs f(K/Nt ohs )} (0.28) 

where the function f(k) resembles a free energy density: it gives the probability of 
observing a substantial deviation between the measured K and its average (K) . (The 



variance of K was also considered in Rcf. (Merolle et al., 2005). In general, this quan- 



tity contains different information to four-point functions such as Xi{k,^) although 
x(tobs) and X4(k,t) may sometimes be related through scaling argument^.) In some 
kinetically constrained models ( |Merolle et al., 2005[ ), the distribution P{K) has a char- 
acteristic shape, skewed towards small activity, with an apparently exponential tail, 
as shown in the central panel of Fig. 10.101 Further, on estimating f(k) from this plot, 
there is a range of K over which f{k) is non-convex (that is, f"(K) < 0). The be- 
haviour of f(k) away from its minimum describes the properties of rare trajectories in 
the system and their relevance for the liquid behaviour is not clear a priori. However, 



6 In the notation of Eq. 110.71 1 , one may write K = ^ . . Cj (t j , t j + At) so that the variance of K 
contains terms like gij mn = (cj(tj , tj + At)c m (t n , t n + At)). Assuming that gij mn depends on scaling 
variables such as \n — rj\/^4 and (tj — t n )/T a , one may connect the variance of K to the dynamical 
length scale §4 and time scale T a . Such connections are are analogous to the relation l|0.14| l between 
X4(i) and ^4(t), but there is considerable freedom in the scaling ansatz for gijmn, which seems to 
prevent a more direct connection between the variance of K and X4 (t ) • 
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the key motivation of this study was to formulate a thermodynamic approach to the 
statistical properties of trajectories of glassy systems ( |Garrahan and Chandler, 2002[ ). 
Within such a framework, non-convexity of f(K) has a direct interpretation as a 'dy- 
namical phase transition' in the system. 

This interpretation is most easily seen by taking the Legendre transform of f(K) 
to obtain a new 'dynamical free energy' 

iP(s) = - min[.sfc + /(&)], (0.29) 

which describes the response of the system to a field s that biases the system towards 
trajectories with small (or large) activity K . In particular, the effect of the bias is to 
change the average of K from (K) to K(s) = —Nt b s ^ip(s). Then, a non-convex form 
for f(K) results in a jump singularity of K (s) for a specific biasing field s = s* . While 
the field s has no simple physical interpretation, one may view it as a mathematically 
convenient trick for sampling the distribution P(K). 

Turning to the results of this formalism, the key point is that glassy systems 
may exhibit singular responses to the field s, leading to 'ideal glass' states that 
are characterised by values of K that is much smaller than its equilibrium aver- 
age (K). The existence of these phase transitions has been proven in simple mod- 
els ( jGarrahan et ai, 2007|Jack and Garrahan, 2010[ ) and numerical results for Lennard- 
Jones model liquids are also consistent with the existence of such a transition ( |Hedges et at, 2009| ) 
In particular, if the field s is chosen to lie at the putative phase transition point, then 
one may construct the distribution of K in the presence of the field s, P S {K) cx 
P(K)c~ sK . In the presence of a phase transition, P S (K) has two peaks, which corre- 
spond to distinct active (liquid) and inactive (glass) states. An example is shown in 
the right panel of Fig. 10.101 

To summarise then, dynamical heterogeneity is concerned with distributions of 
dynamical quantities, most often through their means and covariances. However, the 
tails of these distributions can reveal information about possible new phases in the 
system, whose structure is very stable and whose relaxation times are very long. This 
leads to the hypothesis that the nature of the dynamically heterogeneous fluid state 
should be interpreted in terms of coexistence between and active liquid and inactive 
'ideal glass' states. 

0.5 Open problems and conclusions 

The aim of this chapter was to review the recent progress in the quantitative analysis 
of dynamical heterogeneities. We showed that the introduction of four-point corre- 
lation functions played an important role, both conceptually and operationally, by 
providing a precise quantitative measure to characterise dynamical heterogeneities. 
These four-point correlations have now been measured or estimated in numerical sim- 
ulations of schematic and realistic models of glass formers, and experimentally on 
molecular glasses, colloids and granular assemblies close to jamming. They have also 
been investigated theoretically within simplified models or within the mode-coupling 
approximation, and have indeed been shown to be critical as the glass transition is 
approached. These four-point correlations are the natural counterpart, for glass or 
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spin-glass transitions, of the standard two-point correlations that diverge close to a 
usual second-order phase transition. These can be considered to be breakthroughs that 
have significantly improved our understanding of the microscopic mechanisms leading 
to glass formation, and that have already spilled over to many different scientific com- 
munities. 

However, it soon became clear that four-point correlation functions are not a 
panacea. It was for example not anticipated that these functions would be delicately 
sensitive to details such as the choice of statistical ensemble or the microscopic dy- 
namics. Second, although these four-point objects give valuable information, they are 
not powerful enough to answer more precise questions about the geometry of the 
structures that carry dynamical heterogeneities, or about the nature of the relaxation 
events (continuous vs. activated). Along the same line of thought, the relation between 
the dynamical correlation length extracted from these four-point points and the intu- 
itive (but not so clearly defined) notion of cooperative relaxation is at this stage quite 
elusive. How many different 'dynamical' length scales does one expect in general? 

We have seen that the study of three-point response functions and non-linear sus- 
ceptibilities allows one to bypass some of the difficulties inherent to four-point func- 
tions. We note that the experimental and numerical situation on that front is much 
less developed, and should be encouraged. The very recent measurement of the non- 



linear dielectric properties of glycerol ( |Crauste-Thibierge et ai, 2010[ ) is a remarkable 
exception. 

Higher-order correlation functions might also contain interesting quantitative in- 
formation about dynamical heterogeneities, but this subject is at this stage totally 
unexplored. It was recently suggested ( Lechenault et ai, 2010[ ) that six-point func- 



tions might provide a way to measure intermittent dynamics and identify activated 
events. Skewness (or kurtosis) might indeed detect that the dynamics is intermittent, 
as one expects if 'activated' events dominate, with a few rare events decorrclating 
the system completely, while most events decorrelate only weakly. More work in that 
direction would certainly be worthwhile. 

Finally, we have not touched upon the problem of out-of-equilibrium dynamical het- 
erogeneities, in particular in the aging regime. This is clearly a very interesting topic, 
for which experimental efforts arc underallocated, although results in this regime might 



be able to discriminate between theories. We refer to (Parisi, 1999 .Castillo et ai, 2003 



Parsaeian and Castillo, 2008|Vollmayer-Lee et ai, 2002|Vollmayer-Lee and Baker, 2006| 



El Masri et ai, 2010 Courtland and Weeks, 2003) and the chapter on aging in this 



book for interesting lines of research on this issue in spin-glasses. 
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